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ABSTRACT 



We report the results of a search for molecular oxygen (O2) toward the Orion 
Bar, a prominent photodissociation region at the southern edge of the HII region 
created by the luminous Trapezium stars. We observed the spectral region around 
the frequency of the O2 Nj = 83 - 12 transition at 487 GHz and the 54 - 84 transition 
at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel 
Space Observatory. Neither line was detected, but the 3cr upper limits established 
here translate to a total line-of-sight O2 column density < 1.5 x 10^^ cm~^ for an 
emitting region whose temperature is between 30 K and 250 K, or < 1 x 10^^ cm"^ 
if the O2 emitting region is primarily at a temperature of < lOOK. Because the Orion 
Bar is oriented nearly edge-on relative to our line of sight, the observed column 
density is enhanced by a factor estimated to be between 4 and 20 relative to the face- 
on value. Our upper limits imply that the face-on O2 column density is less than 
4 X 10^^ cm"^, a value that is below, and possibly well below, model predictions 
for gas with a density of 10"^ - 10^ cm~^ exposed to a far ultraviolet flux 10^ times 
the local value, conditions inferred from previous observations of the Orion Bar. 
The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice 
is greater than 800 K; (2) the total face-on Ay of the Bar is less than required for 
O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a 
small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its 
peak abundance, which is density dependent, corresponds to a sky position different 
from that sampled by our Herschel beams. 

Subject headings: astrochemistry - ISM: abundances - ISM: individual objects 
(Orion) - ISM: molecules - submillimeter: ISM 
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1. INTRODUCTION 

Searches for interstellar O2 have a long history, but their motivation has evolved with time. 
Prior to the late-1990's, efforts to detect O2 were driven largely by a desire to confirm its pre- 
dicted role as a major reservoir of elemental oxygen within dense molecular clouds and as the 
most important gas coolant - after CO - of cold (r<30 K), modestly dense (n(H2)::^ 10^- 



10^ cm-^) gas (cf. [Goldsmith & Langer||1978HNeufeld, Lepp, & Melnick||1995|). The launch 



of the Submillimeter Wave Astronomy Satellite (SWAS) in 1998 and Odin in 2001, and the sub- 
sequent failure of these observatories to detect O2 toward a large number of sources at levels 
of a few percent of the abundances predicted by equilibrium gas-phase chemical models, have 
forced a shift in emphasis to a re-examination of the oxygen chemistry in dense molecular gas. 
Today, interest in O2 no longer lies in its being a significant reservoir of elemental oxygen or in 
its cooling power. Instead, because the abundance of gas-phase O2 is set by a balance of various 
formation, destruction, and depletion processes thought to affect the broader chemistry in dense 
gas - such as gas-phase reactions, grain-surface reactions, thermal sublimation, far-ultraviolet 
(FUV) photodesorption, cosmic-ray desorption, photodissociation, and freeze out - measures of 
O2 have become an important test of our current understanding of the relative effectiveness of 
these processes. 

The capabilities of the Herschel Space Observatory 's Heterodyne Instrument for the Far- 
Infrared (HIFI; |de Graauw etal.||2010] ) have enabled improved searches for O2 through: (1) its 
high sensitivity, including at 487 GHz - the frequency of the Nj = Sa - 12 transition observed 
previously by SWAS and Odin; and, (2) its broad frequency coverage that permits observations of 
additional O2 submillimeter transitions, some of which are expected to exhibit stronger emission 
than the 83 - I2 line under certain physical conditions. The Open Time Key Program "Herschel 
Oxygen Project" (HOP; Co-PFs P. Goldsmith and R. Liseau) is designed to survey Galactic 
sources with the goal to detect O2 or set meaningful limits on its abundance within these re- 
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gions. Because the effectiveness of the processes that determine the O2 column density depends 
upon the gas density, temperature, and incident FUV flux Go (scahng factor in multiples of the 
average Habing local interstellar radiation field; |Habing||1968| ) among other parameters, testing 



these models requires that the HOP observations include a range of source types, such as dense 
quiescent clouds, outflows and shocked gas regions, and FUV-illuminated cloud surfaces (see, 
for example, [Goldsmith et aL|20imLiseau etaL|2012| ). 



In this paper, we report the results of a deep search for O2 emission toward the Orion Bar, 
a well known ionization front located approximately 2' southeast of the Trapezium stars in 
Orion at the interface of the HII region created by these stars and the dense gas associated 
with the surrounding Orion molecular cloud. The Orion Bar lends itself well to the study of 
FUV-illuminated molecular gas for several reasons, including its nearly edge-on geometry, its 
proximity 420 pc; [Menten et al.|2007[ [Hirota et al.|[2007t |Kim etaLl[2008] ), its relatively high 



density (n(H2) > 3 x 10^ cm"^), and the strong {Go ^10^-10^) external FUV field irradiating 
this gas. The Orion Bar, and sources like it, are of particular interest since the dust grains within 
these regions are predicted to be sufficiently warm that the thermal evaporation of O atoms 
from the grain surfaces is enhanced, resulting in a higher fraction of O in the gas phase and the 
increased production of O2 via gas-phase chemical reactions (O + OH ^ O2 + H). Under such 
circumstances, the O2 column density can be more than a factor of 10 greater than within gas 
exposed to lower (i.e.. Go < 500) external FUV fields (cf. |Hollenbach et al.|2009[ ). The inclusion 



of the Orion Bar within the HOP program was intended to test this prediction. 

The observations and data reduction methods are described in §2 below. In §3, we present 
the resultant spectra and the upper limits to the O2 integrated intensity. In §4, we review the 
excitation conditions within the Orion Bar and the derived limits on the line-of- sight O2 column 
density. In §5, we discuss these limits in the context of recent chemical models that trace the O2 
abundance from the FUV-illuminated cloud surface to the deep interior. 
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2. OBSERVATIONS AND DATA REDUCTION 

The Herschel HIFI observations presented here were carried out using the HIFI Band la 
receiver for the Sa - I2 487 GHz observations and the HIFI Band 2b receiver for the 54 - 84 774 
GHz observations. The 487 GHz observations were conducted on operational day (OD) 291 in 
spectral scan dual beam switch (DBS) mode, while the 774 GHz observations were conducted on 
OD 297 in spectral scan DBS mode and on OD 509 in HIFI single point DBS mode. Eight LO- 
settings were used for both the 487 GHz and the 774 GHz spectral scans to enable the spectral 
deconvolution, and the additional eight single point 774 GHz observations were observed also 
using eight different LO settings. The total integration time (on- source + off- source) /c^r each 
polarization was 0.93 hours for the 487 GHz spectral scan, 0.86 hours for the 774 GHz spectral 
scan, and a total of 4.6 hours for the eight single point 774 GHz observations. The fuU-width- 
at-half-maximum (FWHM) beam sizes were 44.7'' at 487 GHz and 28.2'' at 774 GHz. 

The observed position, a =5h 35m 20.6s, 5 = -5^ 25' 14.0" (J2000), is shown in Fig.[T] 
We applied the total observing time allotted to HOP observations of the Orion Bar to a single 
spatial position - versus multiple positions - in order to achieve the lowest radiometric noise 
and, thus, the greatest sensitivity to weak O2 emission. In the absence of prior information 
about the possible O2 spatial distribution, our choice of sky position was guided by the desire 
to place the 487 GHz and 774 GHz beam centers a distance corresponding to approximately 
8 visual magnitudes into the molecular gas measured from the ionization front, in accord with 
model predictions (see §5 for a full discussion). For an H2 density between 5 x 10^ cm"^ and 
5 X 10^ cm"^, applicable to the interclump medium in the Bar, and Go — 10^, this corresponds 
to a projected angular distance of between 2.4" and 24" from the ionization front. As shown in 
Fig. [T] the selected position places the beams in the center of this range, while the beam sizes 
encompass the full range. The sky position parallel to the Orion Bar was selected to coincide 
with the molecular gas, as delineated by the ^^CO J =3-2 emission (see Fig. [T]), and, for 
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future analysis, one of the positions under present study by another Herschel Key Program. 



The data were processed using the standard HIFI pipehne software HIPE version 7.3 (Ott 



2010| ), spurious signals ("spurs") removed, spectra defringed, spectral scans deconvolved, and 



all data finally exported to GILDAS-CLASS format. Further processing was performed only on 
the Wide Band Spectrometer (WBS) spectra (0.5 MHz channel spacing and 1.1 MHz effective 



spectral resolution) using the IRAM GILDAS software package (http://iram.fr/IRAMFR/GILDAS/), 
including first-order baseline removal, averaging of the 774 GHz spectral scans and frequency- 
aligned single point observations, averaging of the H- and V-polarization spectra, and production 
of separate averages for both frequencies and both sidebands. The frequencies for the line iden- 
tification were extracted from the JPL and CDMS databases ( [Pickett etal.||1998t [Mtiller et"aL 



2005) as well as Drouin et al. (2010) in the case of O2. 



3. RESULTS 

A summary of the identified lines in the HIFI Band la and Band 2b spectra along with 
the observing modes, integration times, and Gaussian fit parameters is provided in Table 1. 
The summed H-hV polarization spectra observed in Band la are shown in Fig. [2} while those 
observed in Band 2b are shown in Fig. [3} With the exception of the H2CI+ chloronium 485 GHz 



spectrum, which is a blend of three hyperfine components (cf. Lis etal.| 2010; Neufeld etal. 



201 ip , all of the detected lines appear well fit by single Gaussian profiles with a common LSR 
line center of 10.68 ± 0.14 km s~^ (la) and individual best-fit FWHM line widths ranging from 
about 1.8 km s"^ to 2.5 km s"^. 

The upper limit to the integrated intensity of the O2 33-I2 and 54-84 transitions is de- 
rived assuming each line is described by a single Gaussian profile, as is the case for the other 
unblended lines we detect toward this position. The rms noise in the O2 83 - 12 487 GHz spec- 
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tmm between LSR velocities of —1 10 km and +25 km - a velocity range within which 
there is no evidence for any spectral features - is 2.62 mK per 0.5 MHz channel. Similarly, 
the rms noise in the O2 54-84 774 GHz spectrum between LSR velocities of —70 km s~^ and 
+30 km s~^ is 2.19 mK per 0.5 MHz channel. The intrinsic O2 line widths along this line 
of sight are unknown; however, we assume they lie between the extremes of 1.8 km s"^ and 
2.5 km s~^ (FWHM) measured for the other unblended lines we detect along this line of sight 
(see Table 1). This leads to 3a upper limits of between 0.0150 and 0.0209 K km s"^ for the 
33 - I2 487 GHz line and between 0.0126 and 0.0175 K km s'^ for the 54 - 34 774 GHz line. 



4. EXCITATION AND LIMITS ON THE O9 COLUMN DENSITY 



The Orion Bar, like many other photodissociation regions (PDRs), displays emission from 
a variety of ionic, atomic, and molecular species best fit by a mix of gas densities and temper- 
atures. The broad picture to emerge is that of a layer consisting of at least two components: 



interclump gas with n(H2) ^3-20 x 10^ cm ^ ( [Hogerheijde, Jansen, & Van Dishoeck||1995 



Wyrowski etal.|[T997t [Simon et ari|1997[ [Marconi etal.||1998l ) surrounding clumps with n(H2) 



10^-10^ cm-^ (Lis & Schilke 2003; Young Owl etal. 2000), which comprise about 10% 



of the mass ( [Jansen etal.||1995| ). Gas temperature estimates similarly vary, depending on the 
species observed and the component giving rise to most of the emission. Within the denser 



well-shielded gas, the gas temperature is thought to range between ^ 50 and 85 K ( [Hogerheijde, 
Jansen, & Van Dishoeck|1995||Gorti & Hollenbach|2002[ ). The gas temperature associated with 



the interclump medium is estimated to be 85 ± 30 K ( [Hogerheijde, Jansen, & Van Dishoeck 
1995[ ), with some gas temperatures associated with the surfaces (Ay ^ 1) of the denser clumps 
ranging as high as 220 K ( jjansen etal.[[1995t [Batrla & Wilson[[2003i [Goicoechea etd7[[20TT] ) . 
There is evidence for an even warmer component (300- 700 K) based on emission from pure 
rotational lines of H2 and far-infrared fine-structure lines of [OI] at 63 and 145 /xm and [CII] 
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at 158 fim ( [Herrmann etal.|[T997| [AUers etaLl|2005| ). This warmer component is believed to 
arise in the gas between the ionization front and the molecular region traced by ^^CO emission 
( [Walmsley etal.|[2000| ). The strength of the FUV field incident on the Orion Bar has been esti- 
mated to be Go — 1 - 4x 10^ based upon the total radiation from the Trapezium stars - and the 
O star 9^ Ori C in particular - the intensity of the far-infrared [C II] and [O I] fine- structure lines 
mapped toward the Orion molecular ridge, the strength of several near-infrared lines whose in- 
tensities have been ascribed to recombinations to highly excited states of CI, and the strength of 
near-infrared NI lines excited by the fluorescence of UV fines ( [Herrmann etal.||T997[ [Marconi 
etaL[[1998i [Walmsley etaL]|2000| ). Given a density of 10^ cm"^ for the bulk of the material 



and a Go of ^ 10^, models predict that the O2 abundance peaks at Ay ^ 8 mag. (cf. Sternberg 
& Dalgarnol|1995[[Hollenbach etal.[[2009l ). At these depths into the cloud, the gas temperature 
is predicted to be 30-40 K ( [HoUenbach etal.[|2009[ ). Thus, in our analysis, we consider the 
possibility that the O2 emission could arise in gas with temperatures anywhere between 30 K 
and 250 K. 

The weak line flux of the O2 magnetic dipole transitions makes them highly likely to be 
optically thin. Under the assumption that the O2 emission uniformly fills the HIFI beam, the 
observed integrated intensity in a given transition is: 



Tmh dv 



hc^ 



A. N{0,) U = 5.15 X 10-^ ^i^^^ (K km s"^) , 



(1) 



^GHz 



where T^b is the main beam temperature, v is the line frequency (and vqyiz is the line frequency 
in GHz), A^i is the spontaneous decay rate between the transition upper level, u, and lower 
level, I, N(02) is the total O2 column density in cm"^, and fu is the fractional population in the 
transition upper level. The conversion between main beam and antenna temperature makes use 



of the efficiencies reported in Roelfsema etal. (2012) 
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To determine the fractional population of the transition upper state, fu, the excitation of 
the lowest 36 levels of O2, corresponding to a maximum upper-level temperature of 1141 K, 
was computed under the large velocity gradient (LVG) approximation. The spontaneous decay 
rates are those of [Drouin et al.| ( [20T0| ) and the coUisional rate coefficients are those calculated 
by |Lique| ( |2010| ) for He-02 collisions, multiplied by 1.37 to account for the different reduced 
mass when H2 is the collision partner. For molecular hydrogen densities > 3 x 10^ cm"^, both 
the 33 - 12 and 54 - 34 transitions are close to (or in) LTE and the values of fu depend essentially 
only on the temperature. Fig. |4] shows the resulting contours of integrated antenna temperature 
for the 33 - I2 transition as functions of the total O2 column density and gas temperature between 
30 and 250 K. Similarly, Fig. [5] shows the corresponding results for the 54-34 transition. 

Of the two O2 lines searched for here, an examination of Figs. |4] and [5] shows that our 
measured upper limits to the 54 - 34 774 GHz integrated intensity place a more stringent limit 
on the maximum O2 column density for Tgas > 35 K (and comparable limits to that set by the 
487 GHz fine at Tgas ^ 30 K). Specifically, assuming the emission fills the beam, the total line- 
of-sight O2 column density must be less than 1.5 x 10^^ cm"^ (3(7). If the O2 abundance peaks 
within the cooler well- shielded gas, for which Tgas - 100 K, the upper limit to the total O2 
column density is less than 1 x 10^^ cm~^ (3cr). 



5. DISCUSSION 

O2 is produced primarily through the gas-phase reaction O + OH ^ O2 + H and is destroyed 
by photodissociation for the cloud depths of interest here. Thus, the O2 abundance is expected to 
peak where the FUV field has been heavily attenuated and where both the gas-phase O and OH 
abundances are high which, in externally FUV-illuminated clouds, is predicted to occur within 
a relatively narrow (i.e., a few Ay deep) zone centered at an Ay ^ 9 mag. from the cloud surface 
(cf. |Hollenbach etal.||2009| ). The proximity of this zone to the surface and the range of depths 
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over which the peak abundance occurs are governed by several important processes. Near the 
cloud surface, where the FUV field is largely unattenuated, the equilibrium O2 abundance is 
low owing to a high photodissociation rate. Beyond a few Ay into the cloud, the FUV field is 
attenuated, the photodissociation rate reduced, and a region of peak O2 (and H2O) abundance is 
attained. 

Within most clouds with Go < 500, the path to O2 formation is believed to start with the 
formation of water ice, H20ice, on grains, which occurs when O atoms strike and stick to grains 
long enough to combine with an accreted H atom to form OHice and then H20ice. Within this 
region the FUV field remains strong enough to photodesorb H2O from the ice mantles and sub- 
sequently photodissociate these molecules, creating sufficient gas-phase O and OH to produce 
O2 by the gas-phase chemical reaction above. Deeper into the cloud (i.e., greater Ay), the FUV 
field is almost completely attenuated and the gas-phase OH and H2O produced through the pho- 
todesorption and photodissociation of H20ice drops significantly; most O atoms that then strike 
dust grains and form H20ice remain locked in ice as long as the grain temperature is < 100 K. 
Over time (^ 10^ years), this process greatly reduces the gas-phase atomic oxygen abundance 



and suppresses the formation and abundance of O2. Hence, in the model of [HoUenbach etal. 
( |2009| ), the steady-state abundance profile of O2 (and H2O) resembles an elevated plateau that 
peaks at an Ay < 6 for gas with n(H2) = 10^ - 10^ cm"^ and Go < 500. 

For regions subject to a Go greater than ^ 500, such as the Orion Bar, the scenario above 
is altered and, for several reasons, the peak O2 abundance is higher and occurs at a higher Ay. 
First, the high FUV field absorbed at the cloud surface leads to a high infrared field that keeps 
the grains warm, even deep within the cloud. For Go = 10^, Tgr 40 K to Ay ^ 8, resulting 
in a significant fraction of the O atoms being thermally desorbed from the grains before they 
can form H20ice and leading to an increase in O in the gas phase. Second, the higher grain 
temperature also reduces the freezeout of such oxygen-bearing species as OH and O2, further 
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increasing the amount of elemental O in the gas phase. Finally, the attenuated FUV flux at the 
higher values of Ay lowers the photodestruction rates, allowing O2 to survive to greater cloud 
depths. The combined result of these effects is a peak O2 abundance about 3 times higher, and 
a total O2 column density more than 10 times greater than for comparably dense gas exposed 
to Go ^ 500. This result is reflected in the detailed calculations presented in [HoUenbach et al. 



( |2009[ ) and shown in Fig. |6} which is adapted from their paper. For this reason, the Orion Bar 



was considered a promising source for our attempts to detect O2 emission. 

From Fig. [6} it would appear that the upper limits on the total O2 column density established 
here are not in serious disagreement with the model predictions. However, the results shown 
in Fig. |6] apply to a gas column perpendicular to the face of a planar cloud. This is not the 
geometry of the Orion Bar, which has often been described as an edge-on PDR, though its true 
structure has been the subject of some study and debate. For example, based on millimeter and 



submillimeter line observations, Hogerheijde, Jansen, & Van Dishoeck (1995) and Jansen et al. 



(1995) propose a model in which the Bar has a tilt angle, a, of ^ 3° from edge-on, resulting 



in an increase in the line-of- sight column density (beyond what would be measured for a face- 



on geometry) by a factor (sina)"^, or almost 20. Alternately, [Walmsley et al.| ( |2000| ) find that a 
cylindrical model, in which the axis is in the plane of sky and the radius is 0.3pc, best reproduces 
the observed spatial distribution of the fluorescent OI 1.317 /xm emission. In this scenario, the 
average geometrical enhancement of the line-of-sight depth into the Bar versus the face-on depth 
is about 5. Finally, [Neufeld et aL] ( |2006| ) find that a geometrical enhancement factor of ^4 is 



required to reconcile observed and predicted C+ column densities. 

The 3(7 upper limit to tht face-on O2 column density can thus be inferred from our line-of- 
sight values to be 1.5 x 10^^ sin a cm"^, or 1.0 x 10^^ sin a cm"^ for Tgas ^ 100 K. (We note that 
these upper limits are derived assuming the intrinsic O2 FWHM line width is 2.5 km s~^; if the 
intrinsic width is closer to the lower end of the observed range, i.e., 1.8 km s~^, the face-on O2 
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column density upper limits are further reduced by a factor of 1 .4.) For gas densities < 10^ cm~^, 
which applies to most of the gas in the Bar, this is to be compared with a total predicted face-on 
O2 column density of > 7 x 10^^ cm~^, as shown in Fig.[6| with most of this column occurring 
inside a layer of peak O2 abundance with a width corresponding to approximately 2 magnitudes 
(see Fig. [7]), or a linear size of ^ 1.9 x lO^^/n^ cm, where = n(H2)/[10^ cm"^]. Viewed from 
a distance of 420 pc, this zone of peak O2 emission would subtend 3[(l/n^ + 162.4^ sin a]^^ 
where £ is the physical length of the Bar in parsecs. For £ 0.6 pc (cf. [Jansen etaLl|1995| ) and 



77.5 c^l, a>6^ would result in O2 emission that fills the Herschel/HIFI beam at 774 GHz, though 
a minimum geometric enhancement factor of 4, derived from other observations, suggests that 
a does not exceed 15^. However, these tilt angles imply an upper limit to the face-on O2 column 
density between 1.6 x 10^^ cm"^ and 3.9 x 10^^ cm~^, which is below, and in some cases, 
significantly below that predicted by theory. 

For £ 0.6 pc and I, but a < 6^, the O2 layer no longer fills the 774 GHz beam. 
Although the peak O2 column density within the beam will continue to increase for angles less 
than 6^, the beam filling factor will decrease. These two effects offset exactly, and the beam- 
averaged O2 column density will remain the same for all tilt angles less than about 6^. Since the 
O2 emission is optically thin, the line emission will likewise remain constant within the under- 
filled beam. In this case, the geometrical enhancement factor would be ^ 10, and the upper limit 
to the face-on O2 column density remains below that predicted. Therefore, we conclude that 
Bar geometry cannot account for the discrepancy between theory and observations. 

What, then, can account for the discrepancy? The amount of O2 produced in externally 
FUV-illuminated dense gas depends on several factors, which we examine below: 

Thermal evaporation: As noted earlier, the dwell time of an O atom on a grain surface can have 
a considerable effect on the O2 abundance, particularly when this time becomes less than the 
time to combine with an H atom on the surface. The timescale for thermal evaporation of an 
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O atom is approximately 9x10 exp [800 K / Tgr] seconds, where 800 K is the adsorption 
energy of O to water-ice ( [Hasegawa & Herbst||1993| ) that applies to van der Waals binding to a 



chemically saturated surface. It is possible that the binding energy is greater than 800 K, which 
would increase the grain temperature, and thus the Go, required to thermally desorb O atoms 
on short timescales and produce the jump in the total O2 column density for Go Z 500 seen in 
Fig. [6} If, for example, the O adsorption energy was 1600 K, grains as warm as ^ 42 K - the 
expected dust temperature at high Ay in a Go — 10^ field - would, on average, retain their O 
atoms long enough to form H20ice, thus delaying the Go > 500 rise in O2 column density seen 
inFig.|6|until Go> 10\ 

Photodesorption yield ofH20from a grain surface, YH2O' The abundance (and column density) 
of O2 depends on the gas-phase abundance of O and OH, the latter being produced primarily 
through the photodissociation of H2O, much of which is either photodesorbed from grains or 
produced via the dissociative recombination of gas-phase H3O+. At high Go (and Tgr > 20 K), 
short O-atom dwell times on grains suppress the formation of OHice and H20ice. However, 
even though it is not formed on the grain surface in a high-Go environment, H2O formed in the 
gas phase via H3O+ dissociative recombination will be depleted through freezeout onto grains 
and will remain locked in H20ice for as long as rgr< 100 K. Since the quantity of OH and H2O 
returned to the gas phase as a consequence of H20ice photodesorption scales with Yh^o^ the total 
O2 column density likewise scales with Yh^o^ as is seen in Fig. [6j A value for Yh^o less than 
10"^ would help to reconcile theory and observation. However, fits to the SWAS and Odin H2O 
data ( [HoUenbach etal.|[2009| ) as well as theoretical simulations and laboratory measurements 



dAndersson & van Dishoec"k||2008i |Arasa etal.||201li [Westley etal.1|l995a|bi |6berg etal.||2009| ) 



suggest, if anything, that the appropriate value of Yh^o is greater than 10 ^. 

Grain cross-sectional area (perH): The equilibrium O2 abundance in the Ay range of maximum 
O2 abundance scales as (Yh2o)^o-h, where an is the grain cross-sectional area per H nucleus. 
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Therefore, lowering an will decrease the O2 column density, bringing model and observation 



into closer agreement. For an "MRN" (Mathis et al. 1977 ) grain size distribution rigr (a) oc a 



where a is the grain radius, an ^ 2 x 10 cm^ for an assumed gas-to-dust mass ratio of 100 
with grains ranging in radii between a minimum, amin, of 20 A and a maximum, amax, of 2500 A 
(the standard value in HoUenbach etal.|[2009 ). Grains with amin < 20 A will be cleared of ice 



mantles by single photon heating or cosmic rays and, thus, are not significant ice reservoirs. 
Because an oc (amin • <^max) in order to lower the value of an while preserving the total 
mass in grains, either or both amin and a^ax must increase, such as through coagulation. For 
example, a reduction in an, and thus the face-on O2 column density, by at least a factor of 2 
could be achieved if the minimum grain radius were to increase to > 80 A. 

Alternately, the buildup of an ice mantle, which can increase the radius of grains by as much 
as ^ 50 A, will increase the value of an- For values of Go of ^ 10^ applicable to the Orion Bar, 
grain temperatures are expected to be ^ 40 K, which is high enough to inhibit ice formation via 
surface reactions (absent a higher O adsorption energy); however, water formed in the gas phase 
via the reaction H3O+ + e" ^ H2O + H can still freeze out and form an ice mantle. Toward 
Orion, there is evidence for a departure from the assumed gas-to-dust mass ratio of 100, which is 
consistent with the buildup of ice mantles (see, for example, [Goldsmith, Bergin, & Lis|1997| ). In 
addition, there is evidence for a deficiency in small grains and for grain growth, possibly due to 



radiation pressure, the preferential evaporation of small grains, and coagulation (e.g., Cesarsky 



etal.|2000HPellegrini etal.|2009HShaw etal.|2009D . The net effect of lowering an through these 



processes, and increasing an through the accumulation of an ice mantle, is unclear in a high-Go 
environment like the Orion Bar. 

Beam position: For an interclump H2 density between 5x10^ cm"^ and 5x10^ cm"^ and Go = 
10^, the peak O2 abundance is predicted to occur at a face-on depth into the cloud corresponding 
to an Ay ^ 8 (see Fig. [7]). Thus, the linear distance from the Ay = surface, which we assume 
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is the prominent ionization front, to the depth of peak O2 abundance is ^ 7.6 x 10^^/n(H2) cm. 
For an assumed distance of 420 pc, the angular separation between the ionization front and the 
position of peak O2 abundance (and column density) is then 1.5 Av/[n(H2)/10^] arcseconds, 
where Ay is the face-on depth of the O2 peak abundance in magnitudes. Thus, an interclump 
H2 density of 10^ cm~^ should produce O2 emission that peaks ^ 12^' from the ionization front 
and close to the center of the observed sky positions (see Fig. [T]). However, if the interclump 
density is more than a factor of 2 different from 10^ cm"^ - values that remain within the range 
of density estimates for the interclump medium - then the peak O2 abundance is predicted to 
fall to either side of the observed beam center position. 

Finally, we note that the inferred peak line-of-sight H2 column density, A^(H2), applicable to 



the interclump medium toward the Orion Bar is estimated to be 6.5 x 10^^ cm~^ ( Hogerheijde, 
Jansen, & Van Dishoeck||1995| ). If the geometrical enhancement factor is > 10, as would be the 
case for a tilt angle < 5.5^, this would imply a face-on H2 column density of < 6.5 x 10^^ cm"^, 
corresponding to a total Ay through the Bar of about 7. If the face-on extinction through the 
Orion Bar is indeed this low, then the attenuation of the Go ^ 10^ field is not sufficient to allow 
O2 to reach its peak abundance and the total O2 column density will be less than predicted by 
HoUenbach et al.| ( [2009| ), whose total column densities are based upon cloud depths correspond- 



ing to Ay > 10. This is illustrated in Fig. [7} which shows both the profile of O2 abundance versus 
Ay and the cumulative O2 column density to a given Ay, computed using the model described in 



HoUenbach et al.| ( |2009| ) for the conditions appropriate to the Bar interclump medium. At a depth 
corresponding to an Ay of 7, the predicted face-on O2 column density remains < 3 x 10^^ cm"^, 
well below the limits set here. 

The clumps known to exist within the Bar do possess higher H2 densities (i.e., 10^ - 10'' cm~^) 
and column densities (i.e., > 10^^ cm~^; |Lis & Schilke||2003 ) and would provide the necessary 



FUV shielding to allow O2 to reach its full predicted abundance. Such conditions help to recon- 
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cile observation and theory in two ways. First, as shown in Fig. [6} the predicted total O2 column 
densities decrease with higher H2 densities. Thus, the total O2 column density is predicted to 
be lower if the O2 emission arises primarily from within the dense clumps rather than the sur- 
rounding lower density interclump medium. Second, interferometric observations indicate that 
the dense clumps within the Bar typically subtend angles of between A" and 8^' (see, for exam- 
ple, |Lis^^chiIte|2^ and thus provide a natural explanation for why the beam filling factor 
of O2 emission could be less than unity. However, whether the correct explanation for what we 
observe is that O2 emission originates preferentially within the dense clumps, and is suppressed 
within the Ay ^ 7 interclump medium, and with both gas components governed by the processes 
described in |Hollenbach et al.| ( |2009| ), will depend on how well this model reproduces the wealth 
of new lines being detected toward the Orion Bar by Herschel. 



6. SUMMARY 

1. We have conducted a search for O2 toward the Orion Bar, carrying out deep integrations 
around the frequencies of the Nj = 33-I2 and 54-34 transitions at 487 GHz and 774 GHz, 
respectively. Neither line was detected, but sufficiently sensitive limits on their integrated inten- 
sities were obtained to test current models of molecular gas exposed to high fluxes of FUV radia- 
tion - i.e.. Go ^ 10^. In particular, we infer a total face-on O2 column density of <4 x 10^^ cm"^, 
assuming a Bar geometry in which the line-of-sight depth is more than 4 times greater than its 
face-on dimension. This column density is at least 2 times less than that predicted by the model 



of [HoUenbach et al.] ( |2009| ) for the densities, temperatures, and Go appropriate to the Orion Bar. 

2. The discrepancy between the model predictions and our observations would be reduced, 
if not eliminated, if the adsorption energy of atomic oxygen to wate-ice were greater than 800 K, 
and possibly as high as 1600 K. A lower value for the photodesorption yield for H2O would 
help, but is not supported by fits to other astronomical data or recent theoretical calculations 
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and laboratory measurements. A lower grain cross-sectional area per H, such as might occur 
through grain coagulation, radiation pressure, or the preferential destruction of small grains, 
would lower the O2 column density, but it is unclear whether these grain properties apply within 
the Orion Bar. 

3. If the total face-on depth of the interclump medium within the Orion Bar corresponds 
to an Ay < 7, then photodissociation will reduce the O2 column density to values below our 
detection limit. Clumps embedded within the Bar would offer sufficient shielding to enable the 
buildup of higher O2 abundances and column densities in accord with model predictions, while 
the small filling factor of these clumps would reduce the O2 line flux to levels consistent with 
our upper limits. 

4. If the total face-on depth of the interclump medium within the Orion Bar corresponds to 
an Av > 8, it remains possible that most of the O2 emission may have been missed. In particular, 
since the gas density affects the angular separation between the ionization front and the face-on 
depth into the Bar at which the O2 abundance is predicted to peak, interclump H2 densities much 
different than the assumed value of 10^ cm~^ could result in the position of peak O2 abundance 
and column density occurring to either the northwest or southeast of the position we selected. 

Only further modeling, including predictions for other species, can establish which, if any, 
of the above possibilities is most likely to resolve the present puzzle. 

HIFI has been designed and built by a consortium of institutes and university departments 
from across Europe, Canada and the United States under the leadership of SRON Netherlands 
Institute for Space Research, Groningen, The Netherlands, and with major contributions from 
Germany, France and the US. Consortium members are: Canada: CSA, U. Waterloo; France: 
CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, MPS; Ireland, NUI Maynooth; 
Italy: ASI, IFSI-INAF, Osservatorio Astrofisico di Arcetri-INAF; Netherlands: SRON, TUD; 
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Poland: CAMK, CBK; Spain: Observatorio Astronomico Nacional (IGN), Centro de Astrobi- 
ologa (CSIC-INTA). Sweden: Chalmers University of Technology - MC2, RSS & CARD; On- 
sala Space Observatory; Swedish National Space Board, Stockholm University - Stockholm Ob- 
servatory; Switzerland: ETH Zurich, FHNW; USA: Caltech, JPL, NHSC. We also acknowledge 
the effort that went into making critical spectroscopic data available through the Jet Propulsion 



Laboratory Molecular Spectroscopy Data Base ( ,http://spec.jpl.nasa.gov/| ), the Cologne Database 



for Molecular Spectroscopy (http:// www.astro.u ni-koeln.d e/cdms/| and [Muller et aL]|2005| ) and 
the Leiden Atomic and Molecular Database (|http://www.strw.leidenuniv.nl/^moldata/|and [Schoi| er 
et al.|2005] ). Finally, it is a pleasure to acknowledge useful discussions with Dr. Edwin Bergin. 
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Fig. 1 . Position of the HIFI 44.7'' and 28.2'' beams at 487 GHz and 774 GHz, respectively, super- 
posed on an HST image of the Orion nebula ( |0'Dell & Wong||1996[ ). Also shown are contours 
of ^^C O J = 3 - 2 integrated intensity for a portion of a larger map obtained by Lis & Schilke 
( 2003 ), with intensities in K km s"^ noted. The HIFI beams are centered at a = 05^35"^20?6, 
6 = -05^25' 14" (J2000), toward the surface layers of the FUV-illuminated Orion Bar where 
the O2 emission is predicted to peak. 
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Fig. 2. Averaged H and V polarization spectra obtained in HIFI Band la toward the Orion 
Bar, ordered by rest frequency, with the Gaussian fits given in Table 1 superposed in red. Also 
indicated is the energy of the upper level for each transition, in Kelvins. The H2CI+ line is 
a blend of three, partially resolved, hyperfine components (see Table 1). An LSR velocity of 
10.7 km s"^ is denoted with a vertical dashed line. 
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Fig. 3. Averaged H and V polarization spectra obtained in HIFI Band 2b toward the Orion 
Bar, ordered by rest frequency, with the Gaussian fits given in Table 1 superposed in red. Also 
indicated is the energy of the upper level for each transition, in Kelvins. An LSR velocity 
of 10.7 km s~^ is denoted with a vertical dashed line. The frequency of the ^^CO J =7-6 
transition was near the edge of the band and, thus, the spectrum does not cover the full LSR 
velocity range of the other lines. The O2 spectrum has been truncated to remove features in the 
other sideband. 



Page 27 




5 10 15 20 25 

O2 Column Density (cm"2)/l0^5 



Fig. 4. Contours of integrated antenna temperature, in K km s~^, under a Gaussian line profile 
versus line-of-sight O2 column density and gas temperature assuming an aperture efficiency of 
0.7. The shaded area bounds the range of upper limits to the O2 487 GHz integrated intensity 
assuming the intrinsic O2 line FWHM is 1.8 km s"^ (0.0150 K km s"^) or 2.5 km s"^ (0.0209 
Kkms-i). 
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Fig. 5. Contours of integrated antenna temperature, in K km s~^, under a Gaussian line profile 
versus line-of-sight O2 column density and gas temperature assuming an aperture efficiency of 
0.7. The shaded area bounds the range of upper limits to the O2 774 GHz integrated intensity 
assuming the intrinsic O2 line FWHM is 1.8 km s"^ (0.0126 K km s"^) or 2.5 km s"^ (0.0175 
Kkms-i). 
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Fig. 6. Predicted total O2 column density perpendicular to the ionization front as a function of 
Go and n(H+2H2) for H2O photodesorption yields of 10~^ (solid lines) and 10"^ (dashed line). 
The results shown assume a cloud thickness sufficient to encompass the zone of peak abundance 
(after [HoUenbach etal.||2009|). The range of Go that applies to the Orion Bar is shown in the 



shaded region. The horizontal dotted line denotes the upper limit to the O2 column density 
established here, i.e., 1.5 x 10^^ cm"^, divided by a geometrical enhancement factor of 4. 
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Fig. 7. Top panel: Abundance of O2 as a function of face-on depth into a cloud, measured in Ay, 
for a cloud with nn = n(H+2H2)= 10^ cm-^ and 10^ cm-^ exposed to a FUV field of Go = lOl 
This result was computed using the model described in poUenbach et al.] ( |2009| ) assuming their 
"standard" model parameters, except for those noted here. An H2O photodesorption yield of 
10~^ is assumed. The gas and dust temperatures throughout the cloud are calculated self-con- 
sistently in the HoUenbach et al. code, which predicts a gas temperature of 33 K, and a dust 
temperature of 42 K, at the depth of the peak O2 abundance above. Bottom panel: Cumulative 
face-on column density of O2 integrated from the cloud surface to a given depth, in Ay, for the 
abundance profile shown in the top panel. 



